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Generalized Einstein Aggregation Operators Based on the Interval 
Neutrosophic Numbers and Their Application to Multi-attribute 

Group Decision Making 



Abstract: Based on the Einstein operator, the operational rules of interval neutrosophic numbers are 
defined, according to the combination of Einstein operations and generalized aggregation operators, 
the interval neutrosophic generalized weighted Einstein average (INGWEA) operator, interval 
neutrosophic generalized ordered weighted Einstein average (INGOWEA) operator and interval 
neutrosophic generalized hybrid weighted Einstein average (INGHWEA) operator are proposed . And 
some desirable properties of these operators including idempotency, monotonicity, commutativity and 
boundedness are studied, and some special cases including the interval neutrosophic numbers weighted 
Einstein averaging operators and the interval neutrosophic numbers weighted Einstein geometric 
operators are analyzed in the meantime. Furthermore, to address the multiple attribute group decision 
making (MAGDM) problems in which attribute values are interval neutrosophic numbers oriented, a 
group decision making method based on generalized Einstein aggregation operators is developed. 

Keywords: interval neutrosophic number; Einstein operator; generalized averaging operator; multiple 
attribute group decision making 



1. Introduction 

Group decision theory and method have comprehensive application in many fields for instance 
political, economic, military, science, culture and so on. But it is easier to use linguistic information to 
evaluate alternatives and solve the group decision making problems for decision makers. Firstly, Zadeh 
[1] 's fuzzy set (FS) theory proposed in 1965 was a useful tool to deal with fuzzy information which 
has acquired rapid development and a wide range of application in fuzzy multi-attribute decision 
making (MADM) problems. However, the shortage of it is that it has no non-membership other than 
membership . Further Atanassov [2] 's intuitionistic fuzzy set (IFS) which added a non-membership 

function extend the FS, that is to say, the intuitionistic fuzzy sets can be expressed by T A (x) and 

F a (x) which means membership (or called truth-membership) and non-membership (or called 

falsity-membership) respectively. It has made contributions to multi-attribute decision making (MADM) 
and multi-attribute group decision making (MAGDM) [3-5]. Atanassov and Gargov [6], Atanassov [7] 
further developed the interval-valued intuitionistic fuzzy set (IVIFS) and defined some basic 
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operational rules, the most important is these rules extended both the truth-membership function and 
falsity- membership function to interval numbers. However, IFSs and IVIFSs, which can only deal with 
incomplete information and inconsistent information environment which exist commonly in real 
occasions. For example, during a voting process, that thirty percent vote “yes”, twenty percent vote 
“no”, ten percent give up and forty percent are undecided. It can not be presented by the IFS, because 

the indeterminacy ( \-T A (x)-F A (x) ) does not included in IFSs. Seeing that this situation, 

Smarandache [8] presented the concept of the neutrosophic set (NS) which is a generalization of IFS 

and have three parts: T A ( X) , / A (x) and F A ( X) which mean truth-membership, 

indeterminacy-membership and falsity-membership respectively. Obviously, NS is useful to represent 
opinions more accurately and detailedly for decision makers. And we should know that the 

indeterminacy is quantified definitely and T A (x) , I A (x) and F A (x) are independent absolutely in NS. 
Such as the case mentioned above, it can be showed as x(0.3, 0.4, 0.2) from the point view of NS. 

So far, a lot of researches on NSs mostly contain the following: The single valued neutrosophic set 
(SVNS) considered as a case of the neutrosophic set proposed by Wang et al. [9] .The interval 

neutrosophic set (INS) in which T A ( x) , I A (x) and F A (x) were extended to interval numbers defined 

by Wang et al. [10]. The correlation coefficient and weighted correlation coefficient of SVNSs 
developed by Ye [11] who also proved that the correlation coefficient in SVNS is the generalization of 
the cosine similarity degree. The similarity measures between INSs in accordance with the Hamming 
and Euclidean distances represented by Ye [12], then he put forward a multi-criteria decision-making 
method which is on the basis of the similarity degree. To this day, what draw people's eyes is the 
characteristics of INS rather than aggregation operators for INS. 

In addition, it is the information aggregation operators that have become an important research 
direction and got considerable research results. Einstein operators which considered as a special type of 
t-norm have been studied and applied to decision making by some researchers. Therefore, on the basis 
of Einstein operations, Wang and Liu [13] developed the intuitionistic fuzzy Einstein aggregation 
operators which the intuitionistic fuzzy Einstein weighted geometric (IFEWG) operator and the 
intuitionistic fuzzy Einstein ordered weighted geometric (IFEOWG) operator included, and have the 
better smooth approximations than the algebraic operators. Zhao and Wei [14] established intuitionistic 
fuzzy MADM methods based on intuitionistic fuzzy Einstein hybrid average (IFEHA) operator and 
intuitionistic fuzzy Einstein hybrid geometric (IFEHG) operator they developed. Guo et al. [15] 
proposed hesitant fuzzy Einstein weighted geometric (HFEWG) operator, hesitant fuzzy Einstein 
ordered weighted geometric (HFEOWG) operator, hesitant fuzzy Einstein hybrid geometric (HFEHG) 
operator, and hesitant fuzzy Einstein induced ordered weighted geometric (HFEIOWG) operator 
extended Einstein operators to hesitant fuzzy sets. Even so, there is nobody starting researches about 




neutrosophic number aggregation operators according to Einstein operations. 
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Seeing that for one thing INS generalizes NS to interval numbers and can indicate indeterminacy 
more close to the actualities for decision makers, for another arithmetic aggregation operators and 
geometric aggregation operators are the special case of generalized aggregation operators. This paper 
tries to represent some generalized Einstein aggregation operators based on interval neutrosophic 
numbers by means of uniting both Einstein operators and generalized aggregation operators which 
include interval neutrosophic generalized weighted Einstein average (INGWEA) operator, interval 
neutrosophic generalized ordered weighted Einstein average (INGOWEA) operator and interval 
neutrosophic generalized hybrid weighted Einstein average (INGHWEA) operator. In the meantime, 
according to generalized parameter we discuss the properties and special cases of arithmetic 
aggregation operators and geometric aggregation operators. Ultimately, we propose a useful method 
for MAGDM problems in which attribute values are interval neutrosophic numbers oriented. 

To accomplish the goal mentioned above, the remainder of this paper is arranged as follows. In 
section 2, some basic concepts such as INS, Einstein operations and generalized aggregation operators 
are briefly reviewed. In section 3, the new operational laws of INSs based on Einstein operators are 
established and the characteristics are discussed. In section 4, interval neutrosophic generalized 
weighted Einstein average (INGWEA) operator, interval neutrosophic generalized ordered weighted 
Einstein average (INGOWEA) operator and interval neutrosophic generalized hybrid weighted Einstein 
average (INGHWEA) operator which use new operational laws aforementioned in section 3 are 
developed, and some desirable properties and special cases are also discussed. In section 5, a decision 
making method on basis of INGHWEA operator for the MAGDM problems in which attribute values 
are interval neutrosophic numbers oriented. In section 6, an example to illustrate the application of 
proposed method is given, and we make a comparison between the developed method and the existing 
methods. Section 7 summarizes the main conclusion of this paper. 

2. Preliminaries 

2.1 The Interval Neutrosophic Set (INS) 

Definition 1 [8], Let X be a universe of discourse, with a generic element in X denoted by x. A 
neutrosophic set A in X is 

A = {.x(r 4 (x),/ A (x),F A (x)) |xe X} (1) 

Where T A is the truth-membership function, I A is the indeterminacy-membership function, and 
F a is the falsity-membership function. T A (x), / A (x) and F A (x) are real standard or nonstandard 
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subsets of JO , 1 + ^ with not necessarily any connection between them. 

There is no restriction on the sum of T A (x),I A (x) and F A ( x) , so 
0-<T a (x) + I a (x) + F a (x)<3 + . 

The neutrosophic set was originated from philosophy, which leaded it to apply to real life 
difficultly. The single valued neutrosophic set (SVNS) as a generalization of classical set, fuzzy set, 
intuitionistic fuzzy set and paraconsistent sets etc. put forward by Wang [9] which was represented as 
follows. 

Definition 2 [9], Let X be a universe of discourse, with a generic element in X denoted by x. A single 
valued neutrosophic set A in X is 

A = MT A (x),I A (x),F A (x))\xe X} (2) 

Where T A is the truth-membership function, I A is the indeterminacy-membership function, and 
F a is the falsity-membership function. For each point x in X, we haveT A (x), I A (x), F A (x) £ [0,1], 
and 0 <T a (x) + I A (x) + F a (x)< 3 . 

As for the real situations, it is difficult to express T A , I A and F A using crisp numbers, and they 

can be expressed by interval numbers. Wang et al. [11] further defined interval neutrosophic sets (INSs) 
shown as follows. 

Definition 3 [11]. Let X be a universe of discourse, with a generic element in X denoted by x. A 
interval neutrosophic set A in X is 

A = {iT A (x),I A (x),F A (x))\xs X} (3) 

where, T A is the truth-membership function, I A is the indeterminacy-membership function, and 
F a is the falsity-membership function. For each point x in X, we ha veT A (x), I A (x), F A (x) Cl [0,1] , 
and 0 < Sup(T A (x)) + Sup(/ A (x)) + Sup(F A (x)) < 3 . 




For convenience, we can use x = ([T L ,T U ],\I L ,I U ],\F L , F U ]) which can be called an 



interval neutrosophic number (INN) to represent an element in INS. 



Definition 4. Let X = ([7f ,7f ],[/f 7^], [if, /•,' ]) and y = ([7f,7f ],[/ 2 ,I 2 ],[F 2 ,F 2 ])be any 
two INNs, then the cosine of included angle between X and y is 



cos(x,y)=- 



Tff +7f?f +l[l 2 +l"l 2 +F l L F 2 L +/f if 



(f 1 ) 2 +(7f 7 ) 2 +( +(tff +(^f f * p 2 L f +W- Ml +W’ 



Obviously, 0 < COs(x, y) < 1 . When y is the ideal solution I=([1,1],[0,0],[0,0]), the bigger 
the cos (x, /) between X and I is, the more consistent the direction between X and I is. In this 
condition. 



COS (x, I) ■ 



7f+7f 

2 (c77 +( t " 7 +( / , t ) 2 +(/," TfT'-T-TT" ) : 



Definition 5. Let X = (\T', f '], [/', [ if, if ]) and y = ([7f,7f ],[/f if], [if , if ])be any 



two INNs, if COS (x, /) < COS (y,I) then x<y . 



2.2 Generalized hybrid weighted aggregation operations 

The generalized weighted average (GWA) operator is a generalization of the weighted average 
operator, which is defined as follows: 



Definition 6 [16]. Let GWA : ( R j — > R , if 



GWA(a { ,a 2 ,---,a n )= ^ WjO. ) 



where W = (w l ,W 2 ,---,W n ) is the weighting vector of the a , (j = 1,2, ■ • • , n) such that 




Wj e [04]. Z Wj — and /l is a parameter such that Ae (— < °°,0) (J (0,+°°) , then we can call 

j=i 

GW A the generalized weighted average operator. 

Definition 7 [15]. Let GOWA : [i? + ) — > R + , if 



GOWA(a 1 ,a 2 ,---,a n )= ^ cop 



where CO = (<y p fit,, • C0 n ) is the weighting vector associated with GOWA , which 

n 

satisfies CO. E [O.l] , I ®; =1 > b J is the j th largest of (i = 1,2, •••,«) and is a 



parameter such that Ae (— 1 °°,0) (J (0,+°°) , then we call GOWA the generalized ordered 
weighted average operator. 



Definition 8 [16]. Let GHWA : ) — > R + , if 



GHWA(a l ,a 2 ,---,a n )= y'coj) 



where CO = ( CO { , C0-, , ■ ■ ■ , (O n ) is the weighting vector associated with GHWA , which 



satisfies CO. E [0.1] . 



b is the j th largest of 



{nwfli) (i = 1 , 2 , • •• ,n) , W = (vVj, W 2 , ■ W n ) is the weighting vector of (a l ,a 2 ,--- ,a n ) such 

n 

that Wj E [0,l],^ =1, and/lis a parameter such that Ae (- 1 °°,0) (J (0,+°°) , then we can 



all GHWA the generalized hybrid weighted average operator. 



2.3 Einstein operations 



In reality, the t -operators are Intersection and Union operators in fuzzy set theory which are 
symbolized by T-norm (T ), and T-conorm ( T ), respectively. Smarandache and Vladareanu [18] 




introduced a generalized union and a generalized intersection of single valued neutrosophic numbers on 
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the basis of a T-norm and T-conorm . 

Definition 9 [17]. Let X — (Tj , /j , Fy) and y = (T 2 , F , F 2 ) be any two single valued neutrosophic 
numbers, then, the generalized intersection and union are defined as follows: 

*n r>r . y = {T(T 1 ,T 2 ),r(I l ,I 2 ),T\F l ,F 2 )) (9) 

xU rr . y = (r\T l ,T 2 ),T(I 1 ,I 2 ),T(F l ,F 2 )) (10) 

where F denotes a T-norm and T a T-conorm. 

Further, on the basis of T-norm and T-conorm, Einstein operations are defined as follows: 



r E (x,y) 



xy 

l + (l-x)(l— y) 



(ID 



. . x + y 

r £ (*,y) = - — - (12) 

1 + xy 

where F denotes a T-norm and T a T-conorm. 

Let X = (Tj , /j , F x ) and y = (T 2 , 1 2 ,F n ) be any two single valued neutrosophic numbers, then, 

the operational rules of single valued neutrosophic numbers based on Einstein operations are as 
follows: 



(i) X ® E y = 



T l+ T 2 



hh 



F,F 2 



.1 + 7J 2 l + (l-/,)(l-/ 2 ) 1 + (1 -Fj)(1-F 2 ) 



(13) 



(ii) x® E y = 



TT 

1 \ 1 2 



/i + / 2 F t + F 2 



i+(i- 7 ;)(i-r 2 ) i + /,/ 2 l +f,f 2 



(14) 



2 J 



(iii) Ax = 



f {\ +r 1 )"-(i-r 1 ) 1 21* 



2 F* 



(1 + 7 ;) +( 1 - 7 ;/ ( 2 -/,) + 1 * ( 2 -F x ) +F* 



,A> 0 



(15) 




1 

2 

3 

4 

5 

6 
7 



f 



(iv) X =\ 



2T, A 



(i+irf-a-iy (i+F^-a-Fy 



{2~tX +T*' {\ + l x y +{\-l x f ' (l + Fy +(l-F 1 ) 



i ’ i 



V 



, 2>0 



(16) 
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3. Einstein operations of interval neutrosophic numbers 

For establishing the Einstein operations of INNs, we should give the following definition in the 
first place. 

Definition 10. Let X = ([7^F^[^/f and y = ([T^,T^l[I^I^l[F 2 L ,F^J) be 

any two INNs, then, the generalized intersection and union are defined as follows: 



X ® r r . >• = ([r(7/ , t 2 ),y{t^ ,t 2 u )1 [r* (7f , i l 2 X r (I* , I v 2 )] , 
[T\F x L ,F^),r\F? ,F"n) 

X © rr . y = ([r* (T, L , T 2 L ), r"(Tl J , T^)i\r(i 1 ,, / 2 L ), ru \’ , / 2 U )], 
[r (F l L ,F 2 L ),r(F 1 u ,f 2 u )]) 



(17) 



(18) 



where T denotes a T-norm and T a T-conorm. 

According to definition 10, Einstein T-norm and T-conorm, the Einstein operational rules of INNs 
are established in the following. 

Let x = ([7^7f ]) and y = ([T^T 2 U ],[I^I^],[F 2 L ,F 2 U ]) be any two 
INNs, then, the operational rules of INNs based on Einstein operations are shown as follows: 



X© £ y 

( 1 ) 



x® E y 

( 2 ) 



rjiL . rji L rjiU . rpU 

1 l ' ^ 2 1 l ' ^ 2 


I L I L 
l \ 1 2 




1 + T'T 2 ’ 1 + T t u T“ 


<N 

“"-I 

1 

1 

+ 

^ (N 

“--H 

1 

>-1 _ 

1 

+ 
r- H 
1 


J7 l J7 l 

1 1 1 2 


f U f U 
1 1 1 2 


\ 


1 + (1 - F/ Kl - F/ ) ’ 1 + (1 - F t u )(\ - F’ J )\ 


) 


t't 2 


TyT 2 


I L + I L I u +I u 

i l ^ i 2 i l ~ J 2 


1 + (1 - T x l )(\ - T 2 ) ’ 1 + (1 - T” )(1 - T 2 ) 


Li+/ 1 l / 2 ’l+zf/^ 


F\ L + F 2 Fy + F 2 " 


\ 




1 + F''F 2 ’ 1 + F' J F 2 


/ 





(19) 



( 20 ) 
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(l + T, l f - (1 - T, L ) A (l + T” ^ - (1 - T \ u ) A 
(l + T , 1 ) 1 + (\-T x L y' (l + T") A +d-T"Y 



(3) 



21 



LX 

1 



(2-lf-f + !{■* ’(2 




2 R 



LX 



(2 -F x l ) + F { la (2 



2f ua "I 




,X>0 



A - 


f 


t^rj-iLX 27^^ ^ 


A — 

(4) 




(2 -T{) X +T™' (l-T?)* +T™ 
(l+irf-a-i^ (1+1? )"-a- 




(i+fii 


tC 

i 

T 




"-(1 -F?f 


M‘) 




(i+f") 


\ A +(i -f^y 



,X>0 



(21) 



( 22 ) 



Theorem 1. Let x = (XT^ ,F^\)&nd 

any two INNs, then 


y = ( \T 2 l , 7f ] , [ I\ , I u 2 ] , [ F 2 L , F 2 U ] ) be 


(1) x® E y-y® E x 


(23) 


(2) x® E y = y® E x 


(24) 


(3) ti(x® E y) = Jlx® E Jiy,Ti> 0 


(25) 


(4) r/ v x® E r/ 2 x = (rj l + r/ 2 )x, tj v tj 2 >0 


(26) 


(5) x’ 1 ' ® E x'* 1 = x r, ' +r>2 , >0 


(27) 


(6) x 1 ® E y n ={x® E y)\ 7j> 0 


(28) 



It is easy to prove the above characteristics in theorem 1, so omitted in here. 




4. Interval neutrosophic generalized weighted Einstein aggregation 
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operators 



Definition 11. LetJt. =([7J L ,7^ / ],[/j',/y],[i^ L ,/^ / ]) (j = 1, 2, ■ • • ,n) be a collection of the INNs, 
and INGWEA : <Y -> Q , if 



ING WEA ( x, , x 2 , • • • , x n ) 



( n Y" 

® E ( w j x j) 



V M 7 



(29) 



where Q. is the set of all INNs, w = (w { , w 2 , • • • , W n ) T is the weighting vector of 

n 

the (Xj , X 2 , • • • , X n ) such that W . e [0,1]. I W j —l, and A is a parameter such that A G (0, +°°) , 

1= t 

then we can call INGWEA the interval neutrosophic generalized weighted Einstein average 
operator. 



Theorem 2. Let X ] = (\T'\T.],\ Ij 



iu i m 






be a collection of the INNs, 



then, the result aggregated from Definition 1 1 is still an INN, and even 



INGWEA(x v X 2 , • ■ -,X n ) = 



VL 



ID II'-" 



VM 



JN J 



m i Ik 






JN L 



m +i n s ; 



yj=i 



T 

1=1 J 



A f 

+ 



Vx\ 






m-m \m +3 m +\m-m 

V l=i M J V J= l M 7 V J= l 1= 1 7 



^ n n yA f n n yA f n n \/X f n n \/A 

if-' i d> -m-m ii ( ' id -m-m 

V 7=1 J =1 J v j= l J= l J V 7=1 7=1 J v 7=1 7=1 

n n \}a f n n \/A f n n \} / Jl f n n \/x 

II* I D> • ID' ID ID I D‘ - m-m 



Vl=l 



1=1 J 



VM 



M 7 



V i= x 



l=i J 






1=1 7 



n^+aflV 

j 


1 


^ n n ^ 

II' ID 

'v m 1=1 7 


/a 


^ n n ^ 

ID' -ID 


l 


^ « n ^ 

ID' ID. 

1 1=1 1=1 J 


Z' 

1 


^ n n ^ 

1 1 •' -ID 

M 7 


Yx 

+ 


^ n n ^ 

II' ID 

W=i l=i 7 


X’' 


I D'' 'ID 

t H 1- 1 / 


+ 


ID' I lv 

l l=i l=i 7 


/a 

Jy 




( 30 ) 
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where = (l-Tf)* +3(T J L ) X , S, L = (2-Tff -(T^) A 

A t „ =(2-T«f +3 (T") a , B Tr =(2 -T") A -{T«) a , Cjf =(l + lrf + 3(1-1*} 

o , = (i ■ +ij) A - a ■ - 1] f , Cy = a + 1* Y + 3 d - ii y , d iV = (i + i u . r-(\- i v . ) 

M f l=(i+F j l ) A +3(1- F j l ) a ,N fL =(i + F j l ) A -( 1-F j l ) a , M FU =(l+Fj , ) A +3(l-Fj , ) A 



N fU = (1 + Fj ) i - (! - Fj ) a . 



Because of the constraints of space, the proof of Theorem 2 is omitted here. 

It is easy to prove that the INGWEA operator has the following three properties. 

(1) Theorem 3 (Monotonicity) 

Let (x i ,X 2 ,---,X n ) and (jp y 2 »'"» J„) be any two collections of the INNs, if X- < y . for 
all j , then INGWEA(x l ,x 2 ,---,x n )< INGWEA( y x , y 2 , ■ ■ ■ , y n ) . 



(2) Theorem 4 (Idempotency) 

Let Xj =X=([T L ,T u l[I L ,I u \,[F L ,F u ]),j = lZ---,n, then INGWEA(x v x 2 ,- ■ ■ , X n ) = X . 

(3) Theorem 5 (Boundedness) 

The INGWEA operator lies between the max and min operators: 
minCxp^, •••,*„) < INGWEA(x l ,x 2 ,---,x ll )< max(x l ,x 2 ,---,xj 

We will discuss some special cases of INGWEA operator according to the value of A as 
follows. 



(1) When A=1 




INGWEA i=1 (x l , x 2 , • • • , x n ) = © p ( WjXj ) 



nw'-nv-w ni'+^f-nd-TT 

7=1 M M El 

n(i+r/)'' + n<i-7-T' n(i+jfr + n<i-i?)- 



w 



2 fi </")■' 



n( 2 -/‘r + ri(/‘)- n( 2 -/“r + n<;»)- 






zn <?,")• 



fi ( 2 - E y ' 1 + fi < F / ) 1 fi ( 2 - E ) ’ + n < E )"' 



The INGWEA operator degenerates into interval neutrosophic Einstein weighted average 
( INEWA ) operator. 



(2) When A-^0 



INGWEA^ 0 (x l , x 2 , • • • , X„ ) = ® F X W j 1 






2Yl(T< J f 



yi^-tjT 1 n(2- y rf +n^)" 



ro^T-na-/^ n( i+ ^r-na-/rr 

7=1 Tfl M Tfl 

n<i + /; )•■> + n a - /‘r ’ ii(i+i" r + no -/? >- 

. 7=1 7=1 7=1 7=1 

ri (i + e )"' -iid- F N' n(i + e p -iia- e >■ 

7=1 M tl M 

ri(i + F/n + ri(>-r‘r rKi+^r+na-f/y 



The INGWEA operator reduces to interval neutrosophic Einstein weighted geometric 



(. INEWG ) operator. 




Definition 12. Let.x^ =([T J L ,Tj r ],[I I j',I^],[F j L ,Fj , J) (j = 1,2, •••,») be a collection of the INNs, 



and INGOWEA : iY -> Q . if 



INGOWEA(x v x 2 ,-,x n )= © (a>j4 U )) 



where 12 is the set of all INNs, C0= ( 6J t , CO-,, - , C0 n ) is the weighting vector associated with 

n 

INGOWEA , such that CO ■ 6 [0, 1 ] , CO- = 1 . (tJ(l ), (J(2), • • • , (J(/7 ) ) is a permutation of 



(1,2, •••,«) which satisfies X a(j]) > X (T(]) for all j , and A is a parameter such that /le(0,+°°), 

then we can call INGOWEA the interval neutrosophic generalized ordered weighted Einstein 
average operator. 



Theorem 6. Let X- = ([72 , Tj ], [Ij , I ], \F- , F. ]) ( j = 1, 2, • • • , n) be a collection of the INNs 
and X > 0 , then, the result aggregated from Definition 12 is still an INN, and even 
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ING O W EA ( x , , x, , • • • , x n ) 




where A pf = (2-7^ +3(7*)* , =(2-lff -(jf) 1 

=(2-Tff +3(jfY , B T „=(2-T' j ’) A -(T^f , C /; -(l + /'-)" +3(1 , 
d i} =(i- +ij) A - (1-ijY , C'^ii+I^Y+Sii-IIY , D iV = <\+i v .f-<\-i v .f , 
M f l = (l+ F/ ) i + 3(1 - Ff) 1 ,N^=(l + Ff ) A - (1 - Ff Y , = a+7^)*+3(l-7f )\ 
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Limited to the space, the proof of Theorem 6 is omitted here. 

Because the input arguments are weighted according to the arguments position in descending 
order for the whole arguments in the INGOWEA operator, so the OJ = (o) r OX, - ■ ■ , OJ n ) is called 

the position weighted vector. Position weighted vector CO = ( G ) x , C0 2 , • • • , CO n ) should be assigned 
according to the real decision-making needs generally. 

It is easy to prove that the INGOWEA operator has the four properties. 

(1) Theorem 7 (Monotonicity) 

Let (Xj , X 2 , • • • , X n ) and ( Tp V 2 > ' " > T„ ) be any two collections of the INNs, if Xj < y . for 
all j , then INGOWEA( x l ,x 2 ,--’,x n )< INGOWEA( y l , y 2 , • • • , y n ) . 

(2) Theorem 8 (Idempotency) 

Let Xj =x = ([T L ,T u ],[I L ,I u l [F l ,F u ]), j = 1, 2, ■ • n , then INGOWEA{x x ,x 2 ,—,x n ) = x. 

(3) Theorem 9 (Boundedness) 

The INGWEA operator lies between the max and min operators: 
minfx, , x 2 , • • • , x„ ) < INGOWEA(x l , x 2 , • • • , x„ ) < max(xj , x 2 , • • • , x„ ) 

(4) TheoremlO (Commutativity) 

Let (Xj,X 2 ,---,X n ) be any permutation of (_y p _y 2 ,---, _y n ) , then 
INGOWEMx, , x 2 , • • • , x„ ) = INGOWEA(y , ,y 2 ,-,y n ). 

Similarly, as for different values of parameter A , some special cases of INGOWEA operator 
will be discussed in the following. 



(1) When A= 1 




n 

INGOWEA A=l (x l , x 2 , ■ ■ ■ , x n ) = ©£ ( 0)jX aU) ) 

7= 1 

[ ri f 1 ■ +y > r - n a ■ - ■ c )*' n h + ■ c r - fi « 

7=1 M M 7=1 

n(i+^ u ,)'' +n<i-^ 0 ,)'' ii(i+c,r +ii<' 

7=1 7=1 7=1 7=1 



rjiU \ ^ 

"OOV 



7=1 



ri(i+y>r +n<'-i'«„)'' ii(i+c„r +ri<i-o'' 

L 7=1 7=1 7=1 7=1 J 

[ 2 n<c/' 2 n<'"o>)"' 1 

= — — — , (35) 

n(2-'«„r +n<o*' n( 2 -c„r 

L 7=1 7=1 7=1 7=1 J 

[ 2 fl(^>)"' AVC,f 1 

M M_ 

ri e-o" + ii<^,)*' riE-o" +n(0"' 

VL 7=1 7=1 7=1 7=1 J 

The INGOWEA operator degenerates into interval neutrosophic Einstein ordered weighted 
average ( IN E 0 WA) operator. 

(2) When /l->0 

n 

INGOWEA^ q (x x ,x 2 , • ••,*„) = ® E x‘°J U) 

7=1 

[ 21W } 

M M 

ri( 2 -^j”'+ri<o*' n( 2 -c.r + ii<c.)“' 

L 7=1 7=1 7=1 7=1 J 

[ n( i+, «y.)'' -iic-o'' n(i+c„r -rto-o** 1 

= ^ ^ , (36) 

n(i+'i<„f +n<i-'« fl )"' n( i+ o"' +n< i -o" 
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The INGOWEA operator reduces to interval neutrosophic Einstein ordered weighted geometric 




( IN E OWG) operator. 
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What we discussed the two interval neutrosophic Einstein operators, including both INGWEA 
operator and INGOWEA operator. One hand the INGWEA operator lays stress on their own 
importance of these arguments, while on the other hand INGOWEA operator emphasizes their 
ordered position importance. However, as for real decision making problems, we need take consider 
the two weighted modes into consideration which indicate different aspects of decision-making 
problems. In order to overcome the weakness that each operator considers only one aspect, 
generalized hybrid averaging operator based on Einstein operations is proposed as follows. 



Definition 13. LetXj =([T j l ,T^\,[Ij,I 1 I\,\F j L ,F^J) (j = 1,2, •••,«) be a collection of the INNs, 



and INGHWEA : Q" -> Q , if 



,1/A 



INGHWEA(x i ,x 2 ,---,x„) = 



© r (^VVj)) 



E 

V 1=1 



(37) 



where Q, is the set of all INNs, CO = (filj , (0 2 , • • • , CO n ) is the weighting vector associated with 

n 

INGHWEA , such that CO e [0, l] , ^ CO j = 1 . w = (\V X , W 2 , • • ■ , W n ) J is the weighting vector of 



j = i 

n 

the(x 1 ,x 2 ,---,x n ) such that Wj G [O.lJ.^W', = 1 . Let .v ; = nw j X ] = (\t'\ /' 7 ], [F'\ F' j ]) , 

l=i 

n is the adjustment factor. (<J(1), (J(2 ), ■••,(7(n)) is a permutation of (1,2, •••,«) which satisfies 

x a(j-i) — x a( D f° r J ’ an d 2- is a parameter such that /le(0,+°°), then we can call 
INGHWEA the interval neutrosophic generalized hybrid weighted Einstein average operator. 



Theorem 11. Let Xj ={[T^T^l[I- (j = 1,2, •••,«) 



be a collection of the 



INNs and X > 0 , then, the result aggregated from Definition 13 is still an INN, and even 
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INGHWEA( x v x 2 ,--,x n ) 




where = (2-7\ L )* +3(7^ , = (2-77 L f -(T^f 

=(2-T”) 1 +1{T?) X , = (2-T") A -{T") x , C /t = (l + 1) 'f + 3(1 - / ) i , 

D,. = (l + / 7 L ) A -(l-/ 7 L / , C /y =(1 + /?)* + 3(1-/;/ , D ;; =(l+/7) A -(l-/7)^ , 
M Ff =(l+F j L ) A +3(\-F j L f , N Ff , =(\ + Ff) A -a-F^ , = Q+Ff) x +3(1-^ , 
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(l +T L ) nWf —(\-T L ) nw ’ (l +T u V i -(l-T u ) n 

N = (1 + F u ) z -(1-F u Y T l - [ 1 ' — T u ~ 1 ’ — 

F * 1 1 ’ J (l+^p+a -T\r r 7 (i+^p+a-^r 



jL _ WP 

J / . _ i \nWi 



t 2 -ijr+w 



„ 2 (F")~' 



i u =- 



2(i, r 



(2 -/yp+ayp 



. F7=- 



2(Fj) 



(2 -Ff)" w ‘ +(Ff) n 



F = 



1 (2 -F”] m> +(F”)" 



Limited to the space, the proof of Theorem 1 1 is omitted here. 



Theorem 12. The INGWEA operator and INGOWEA operator are the special cases of the 
INGHWEA operator. 



(l)Iftf> = 





INGHWEA operator reduces to INGWEA operator. 



(2) If W = 



f i 

V n n 



n 

n) 



then 



INGHWEA operator reduces to INGOWEA 



operator. 



Similarly, when A = 1 , the INGHWEA operator will degenerate into interval neutrosophic 
Einstein hybrid weighted average ( IN EH W A) operator, and when A — > 0 . the INGHWEA 

operator will degenerate into interval neutrosophic Einstein hybrid weighted geometric ( INEHWG ) 
operator. 

According to definition 13, we can know that the INGHWEA operator firstly weights the 
given input arguments, and then reorders the weighted values in descending order and weights these 
ordered arguments. Therefore, the INGHWEA operator can reflect the importance degrees of both 
the given input arguments and their ordered positions. 



5. Multiple attribute group decision making method based on 
INGHWEA operator 



Taking the MAGDM problems based on INNs into account: let A = {A p A 2 ,---, A,,} , 




C — {Cj, C 2 , • ■ •, C (J }be the set of alternatives and attributes respectively, and vv is the weight of the 

n 

attribute C (j = 1,2, , where 0 < U' ; < 1 (j = 1,2, ■••,«) , — 1 . Suppose that 



D={D l ,D 2 ,-,D l } is the set of decision makers, and y t (k = 1,2, ■••,?) is a weight of decision 



maker D, 



withO<y t <1 (k = 1, 2, — , f) , V y k =1. Suppose that X <k> — ( X* 1 ) is the decision 

““ * \ J / rnx.n 



matrix, where xf =([T? k) ,I^ k) ],[F^ k \F^ k) ]) takes the form of the INN and 

y y v y y y ” 

which means that the decision maker D k evaluates the attribute C ; with respect to the alternative^.. 

Then, we can rank the order of the alternatives based on the given information. 

Attributes can be classified into two types including benefit attributes and cost attributes generally, 
and the n attributes may be evaluated in different ways, so what we should do is normalizing the 
decision-making information before applying the INGHWEA operator to the MAGDM problems. 

The method involves the following steps: 

Step 1. Normalize the decision-making information. 

Suppose R (k) — ( E k) ) (k = 1,2, is the normalized matrix of the decision matrix 

XW= ( X r) mx „ (* = U,-,0, where rf = , then the 

normalization method is chosen as follows: 

(1) For benefit attributes: 



where i = l,2,---,m, j = 1, 2, •••,«, k = l, 2 , 



(2) For cost attributes: 






where i = l,2,---,m, j = 1 , 2, •••,«, k = l,2,---,t . 




Step 2. Aggregate the assessment information of each decision maker D k by INGHWEA operator, 
and get the collective information. 
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Step 3. Calculate the comprehensive evaluation value of each alternative. 



r i =([l,T i UlJ i m,F i ]) = INGHWEA(r n ,r i2 ,-,r m ) 



(41) 



(42) 



Step 4. Calculate the cosine of included angle between r (i = 1, 2, • • in) and the ideal solution 

1 = ([ 1 , 1 ], [ 0 , 0 ], [ 0 , 0 ]). 



cos (rj) = 



T, + T, 



^((ZY +{T ,) 2 +Uf +{I ,) 2 +(Z ) 2 +(F i ) 2 ) 



(43) 



Step 5. Rank the alternatives { Aj , A, , • • • , A in } in descending order according to the Definition 5. 

Step 6. End. 

6. An illustrate Example 

Suppose that the example about the air quality evaluation of Guangzhou for the 16th Asian 
Olympic Games (adapted from [18]). The air quality in Guangzhou for the Novembers of 2006, 2007, 
2008 and 2009 were collected as the set of alternatives that is{A p A 2 ,A 3 , A 4 } = [November of 2006, 
November of 2007, November of 2008, November of 2009}. There are three air-quality monitoring 
stations expressed by ( D { , D-, , ) which can be seen as decision makers, and their weight 

y= (0.314,0.355,0.331) . The measured attributes (suppose that their weight vector 

is W = (0.40,0.20,0.40) ) are shown as follows: S0 2 (Cj ) , N0 2 (C 2 ) , PM 10 (C 3 ) . 

The measured values from air-quality monitoring stations under these attributes are shown in 
tables 1, 2 and 3, and they can be expressed by INNs (Note: the original data take the form of interval 
intuitionistic fuzzy numbers, we can get INNs by I L =1 —T u — F u and I u =1 — T L —F L ). 



Table 1 Air quality data from station Z)j 
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C, C 2 c 3 

A ([0.22, 0.31], [0.15,0.55], [0.23, 0.54]) ([0.13, 0.53], [0.1 1.0.67], [0.20, 0.36]) ([0.12, 0.37], [0.07 ,0.48], [0.40, 0.56]) 

A ([0.28, 0.41], [0.10,0.39], [0.33, 0.49]) ([0.33, 0.53], [0.1 1.0.47], [0.20, 0.36]) ([0.12, 0.37], [0.17, 0.58], [0.30, 0.46]) 

A ([0.32, 0.41], [0.15,0.45], [0.23, 0.44]) ([0.43, 0.53], [0.22,0.41], [0.16, 0.25]) ([0.23, 0.45], [0.18, 0.56], [0.21, 0.37]) 

A ([0.39, 0.47], [0.17,0.43], [0.18, 0.36]) ([0.39, 0.53], [0.15,0.34], [0.27, 0.32]) ([0.28, 0.34], [0.43, 0.61], [0.1 1, 0.23]) 

Table 2 Air quality data from station D 1 

c, c 2 c 3 

A ([0.04, 0.21], [0.33,0.61], [0.35, 0.46]) ([0.10, 0.341, [0.21, 0.63], [0.27, 0.45]) ([0.32, 0.37], [0,43,0.55], [0.13, 0.20]) 

A ([0.32, 0.39], [0.22,0.41], [0.27, 0.39]) ([0.03, 0.57], [0.07 ,0.67], [0.30, 0.36]) ([0.16, 0.25], [0.56,0.70], [0.14, 0.19]) 

A ([0.26, 0.37], [0.23,0.53], [0.21, 0.40]) ([0.23, 0.43], [0.42, 0.71], [0.06, 0.15]) ([0.21, 0.35], [0.36,0.68], [0.11, 0.29]) 

A ([0.30, 0.43], [0.22,0.51], [0.19, 0.35]) ([0.28, 0.43], [0.23, 0.41], [0.31, 0.34]) ([0.39, 0.46],[0.37,0.60],[0.01, 0.17]) 

Table 3 Air quality data from station D, 

C, C 2 c 3 

A ([0.25, 0.27], [0.33, 0.52], [0.23, 0.40]) ([0.17, 0.27],[0.33,0.57],[0.26, 0.40]) ([0.21, 0.30], [0.38,0.62], [0.17, 0.32]) 

A 



([0.25, 0.29], [0.32, 0.42], [0.33, 0.39]) ([0.18, 0.46], [0.04,0.39], [0.43, 0.50]) ([0.06, 0.21], [0.49,0.66], [0.28, 0.30]) 
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c, 



c 2 



c 3 



^3 ([0.22, 0.27], [0.42, 0.51], [0.27, 0.31]) ([0.13, 0.37], [0.43, 0.71], [0.16, 0.20]) ([0.11, 0.24], [0.57, 0.75], [0.14, 0.19]) 

A 4 ([0.30, 0.48], [0.07 ,0.61], [0.09, 0.45]) ([0.08, 0.53], [0.23,0.72], [0.20, 0.24]) ([0.32, 0.61], [0.30,0.67], [0.01, 0.09]) 



6.1 The evaluation steps 



The steps are shown as follows: 



(1) Normalize the decision-making information. 



Because the whole measured values are of the same type, so they do not need normalization. 



(2) Aggregate the measured values of each air-quality monitoring station D k by 



INGHWEA operator in formula (41), suppose A = 2 and CO — 



1 1 1 



yb 3 3y 



Then we can get an aggregated matrix R based on three decision matrices measured by 
air-quality monitoring stations, which is shown as follows: 



I[0. 192, 0.264], [0.239, 0.560], [0.267, 0.439]) ([Q 136,0.396], [Q20C1 Q618], [0.243, 0.401]) 
([0.28% 0.369], [Q 196, 0.407], [0.307, 0.417]) ([0215,0.526], [0.067, 0.495], [0.297, 0.399]) 
([0.209, 0.356], [Q2450.497], [0.235, 0.377]) ([Q2930.448], [0.3)5, 0.582], [0.114, QM]) 
( 1 033 1, 0.460|, 1 0. 140 0.5 10|, 1 0. 146 0.383]) ([0284, 0.497], [0201,0.451], [0257, 0.297]) 



([0244,0350], [0237,0546], [0202,0316])^ 
([0.124, 0.283], [0365,0.616], [0.224, 0.289]) 
([0. 192,0.3581, [0.3330.655], [0.1460.271]) 

( 1 0.340 0.494-1, 1 0.300 0.624|, 1 0.022 0. 152|) y 



(3) Calculate the comprehensive evaluation value of each alternative by formula (42), suppose 
A = 2 and CO — 

r x =([0.223, 0.334], [0.234, 0.545], [0.229, 0.372]) r 2 = ([0.234, 0.373], [0.206, 0.490], [0.261, 0.343]) 
r 3 = ([0.252, 0.382], [0.287, 0.553], [0.168, 0.288]) r 4 = ([0.344, 0.494], [0.219, 0.532], [0.077, 0.246]) 

(4) Calculate the cosine of included angle between r 0 = 1,2,- • -,m) and the ideal solution 
I = ([l,l],f0,01,[0,0]) based on formula (43). 



Ill 

3’3’3 




cos* (?;,/) = 0.469, cos*(r 2 ,/) = 0.528, cos*(r 3 ,/) = 0.532, cos* (r 4 , /) = 0.680 
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(5) Rank the alternatives 

According to the comparison method described in Definition 5, the ranking is 
A 4 >- A 3 >- Aj y- Aj , So the best alternative is A 4 , that is, the best air quality in Guangzhou is 
November of 2009 among the Novembers of 2006, 2007, 2008, and 2009. 

6.2 Discussion 

For demonstrating the influence of the parameter A on decision making of this example, we 
change the value of A in the second and third step to rank the alternatives, then get the ranking results 
are shown in Table 4. 

Table 4 Ordering of the alternatives by utilizing the different A in INGHWEA operator 



A 



cos* (?;,/) 



Ranking result 



A = 0.01 



/l = 0.5 



A = l.O 



A = L5 



cos* (Vj , I ) = 0.398 
cos*(r 2 , /) = 0.460 
cos *(r 3 , /) = 0.486 
cos *(r 4 ,/) = 0.637 

cos'Oj, 7) = 0.399 
cos (r 2 ,I) = 0.460 
cos *(r 3 , /) = 0.485 
cos *(r 4 ,/) = 0.636 

cos’Oj,/) = 0.417 
cos (r,,/) = 0.478 
cos *0 3 , 1) = 0.496 
cos *0 4 ,/) = 0.647 

cos’Oj ,1) = 0.443 
cos Ot , 7) = 0.502 
cos 0 3 ,/) - 0.513 
cos*0 4 ,/) = 0.663 



A 4 >- A 3 >- A2 y Aj 



A 4 >- A 3 >- A, ^ Aj 



A 4 Aj A, ^ Aj 



A 4 A 3 >" Aj ^ Aj 
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A = 2.0 



A = 2.3 



A = 2.4 



A = 5 



A = 10 



cos*(/-j, /) = 0.469 
cos *(r 2 ,/) = 0.528 
cos (r 3 , /) = 0.532 
cos*(r 4 ,/) = 0.680 

cos *(rj,7) = 0.484 
cos (r 2 , 1) - 0.5425 
cos (r 3 , /) = 0.5435 
cos *(r 4 , /) = 0.6895 

cos*(r 1 ,7) = 0.489 
cos \r 2 ,I) = 0.5474 
cos*(r 3 , /) = 0.5471 
cos *(r 4 , /) = 0.693 

cos’(r„/) = 0.587 
cos *(r 2 , /) = 0.650 
cos *(r 3 , /) = 0.620 
cos*(r 4 , /) = 0.757 

cos* (Yj , / ) = 0.678 
cos*(r 2 , /) = 0.747 
cos *(r 3 , /) = 0.688 
cos*(r 4 , /) = 0.818 



A 4 >- A 3 >- A, >- Aj 



A 4 >- A 3 >- Aj ^ Aj 



A 4 Aj A 3 A, 



A 4 Aj A 3 A, 



A 4 Aj A 3 >- A, 



As we can see from Table 4, the ordering of the alternatives may be different for the different 
value A in INGHWEA operator. 

(1) When 0 < A < 2.3 , the ordering of the alternatives is A 4 >- A 3 >- A, >- Aj and the best 
alternative is A 4 . 

(2) When A > 2.4 . the ordering of the alternatives is A 4 >- A 1 >- A 3 >- A, and the best 
alternative is A 4 . 

From above analysis, we find that the ordering of the alternatives changes with the value of 
parameter A . To some extent, A can be considered as the mentality of the decision-makers, the more 
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the X is, the more optimistic decision-makers are, and vice versa. Generally speaking, we consider the 
rankings of some special values when X gets 0 and 1 as what we want. We can think the best 

alternative is A 4 in this example. 

For confirming the effective of the proposed method, we utilized the method proposed by Yue [18] 

to solve this illustrate example, and we can get the ranking result A 4 >- A 3 >- Aj >- A, . Because the 

proposed method has a variety of weighted average modes with the change of parameter X , such as the 
geometric average, arithmetic average, quadratic average, and so on. WhenO < X < 2.3 , the ordering 

of the alternatives is A 4 >- A 3 >- Aj >- A x , which is the same as the ranking result by using the method 

proposed by Yue [18]. However, the method proposed in this paper become more flexible because of 
parameter X . 

7. Conclusion 

During the process of solving the real problems, the decision maker always encounters the 
evaluation information of alternatives which is incomplete, indeterminate and inconsistent. Fortunately, 
the interval neutrosophic set (INS) is a better tool to depict this kind of information. In addition, 
Einstein operations meet the typical t-norm and t-conorm and have better smooth approximations than 
the algebraic operators. Therefore, in this paper, considering interval neutrosophic fuzzy information 
and on the basis of Einstein operators and generalized aggregation operators , interval neutrosophic 
generalized weighted Einstein average (INGWEA) operator, interval neutrosophic generalized ordered 
weighted Einstein average (INGOWEA) operator and interval neutrosophic generalized hybrid 
weighted Einstein average (INGHWEA) operator are proposed . Meanwhile, some properties of these 
operators are explored and special cases are discussed with respect to generalized parameter. 
Furthermore, we have given a decision making method based on INGHWEA operator. 

This paper further extends Einstein operators to interval neutrosophic fuzzy information, and 
establishes a group decision-making method based on Einstein operators, which not only enriches 
neutrosophic set theory, but also further expands the using scope of Einstein operators. However, the 
proposed method only solves the MAGDM problems where the attribute values are INNs and the 
weights of attributes and decision makers are real numbers. In further research, we should study the 
MAGDM problems where the weights of attributes and decision makers are INNs, which exist 
commonly in real situations. In addition, it is necessary and meaningful to extend the applications of 
the operator to the other domains such as pattern recognition, fuzzy cluster analysis and uncertain 
programming, etc. 
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